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Abstract. Given the probability measure v over the given region Q C M n , 
we consider the optimal location of a set E composed by n points Q in order 
to minimize the average distance S i-> J n dist (x, S) dv (the classical optimal 
facility location problem). The paper compares two strategies to find optimal 
configurations: the long-term one which consists in placing all n points at once 
in an optimal position, and the short-term one which consists in placing the 
points one by one adding at each step at most one point and preserving the 
configuration built at previous steps. We show that the respective optimiza- 
tion problems exhibit qualitatively different asymptotic behavior as n — * oo, 
although the optimization costs in both cases have the same asymptotic orders 
of vanishing. 



Planning an economic activity is in general an extremely complex problem, where 
a high number of parameters often intervene. In addition, the attitude of the plan- 
ner has to be taken into account: long-term planners take their decision through an 
optimization process over a large time horizon, while short-term planners behave op- 
timizing day-by-day their strategies. Usually, the first kind of behavior is perceived 
as more virtuous and efficient, while the second is seen as easier to implement. 

In this paper we analyze the long-term and short-term strategies in a very simple 
model problem, and we give a way to measure the efficiency of the first versus the 
second. The problem we consider is the so-called location problem which can be 
roughly described as follows. Suppose one has to open a certain given number 
n € N of identical facilities (e.g. plants, shops, distribution centers, cinemas etc.) 
in the given urban region f2 C M. d which will be modeled by a compact convex set. 
If the density of population in Q is given by a known Borel probability measure u, 
then the simplest way to measure how good is the chosen configuration of facilities 
modeled by a set E C ft consisting of at most n points (i.e. #E < n), is clearly to 
calculate the average distance the people have to cover to reach the nearest facility 



Hence the owner of the facilities is interested in locating them in such a way as to 
minimize F. That is why, if he is able to open all the facilities at once, he would 
choose a configuration E = S n solving the following problem. 

Problem 1.1. Minimize the functional E i— » -F(E) subject to the constraints 



We will further refer to such problem as a long-term planning problem since it 
models the optimal choice of facility location so as to satisfy global (long-term) needs 
of the facility owner. If E = E„ solves Problem 11.11 we will denote l n :— F(E n ) 
the respective optimal cost. It is worth mentioning that such a problem has been 
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1. Introduction 
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Sc!l, #E < n. 
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extensively studied (see e.g. [151 El HH] for recent surveys on the subject), but 
nevertheless a lot of interesting and important questions regarding this problem 
still remain open. 

If however the owner of the facilities is unable to open all the facilities immedi- 
ately (say, if he does not have enough financial resources to do that), he will try 
to open them one by one, trying to minimize the average distance functional F at 
each step (i.e. when opening each facility), but, of course, taking into consideration 
the location of facilities already opened at previous steps. Such a short-term opti- 
mization strategy amounts to solving the following problem (in the sequel referred 
to as short-term planning problem). 

Problem 1.2. For each n € N, n > 1, find a set £ = minimizing the functional 
£ i ► F(S) subject to the constraints 

K-i c£c!],#E< n, 

where £q := 0. 

The above model, to the best of our knowledge, can be traced back to [15], where 
the short-term facility allocation strategy is called a myopic allocation policy (see 
also references therein for similar formulations). In the same paper the authors pro- 
pose also several other allocation strategies which can be considered as intermediate 
between the short-term and the long-term ones. 

We will further denote s n := F(E' n ) the optimal cost of the above short-term 
optimization problem. 

In this paper we discuss the asymptotic behavior of solutions to the above two 
problems as n — + oo, by studying the weak limits (in a suitable sense) of optimal 
configurations £„ and Y,' n as well as the optimal costs l n and s n . Namely, we will 
be interested in finding answers to the following questions. 

(A) Find the asymptotic order of l n (resp. s n ) as n — > oo, i.e. find an exponent 
a > such that 

dn- a < l n < C 2 n- a (resp. C x n' a < s n < C 2 n- a ) 

for some positive constants C\ and C 2 and for n sufficiently large. To 
simplify the notation, in the sequel we will write in this case /„ ~ n~ a 
(resp. s n — n~ a ). 

(B) Find precise asymptotic estimates for l n (resp. s n ), i.e. find lim„ n a l n (resp. 
lhn n n a s n ) 1 or just liminf and limsup, should the limit not exist. 

(C) Describe the asymptotic behavior of the minimizers, i.e. find all the weak 
limits, in a suitable sense, of subsequences of minimizers. 

In the pioneering paper [15) some accurate numerical calculations of optimal 
short-term configurations for the case of the uniform density on a line and on a 
two-dimensional square have been provided for rather small values of n. Moreover, 
some deep insights have been formulated regarding the behavior of solutions. Nev- 
ertheless, the above questions have not been explicitly addressed nor even formally 
posed. In this paper we provide rigorous results concerning the nature of the prob- 
lem which partially confirm the insights of [15], and also go further. In Section [2] 
we summarize all the known results on the asymptotic behavior of solutions to the 
long-term optimization Problem 11.11 that we need for the purpose of comparison 
with the short-term optimization Problem 11.21 In Sections [3] and [4] we show that 
although the optimal cost of the short-term optimization Problcm ll.2l has the same 
order of asymptotic expansion as n — ► oo (i.e. the answer to question (A) is the 
same for both problems), nevertheless, the two problems are qualitatively different 
even in the simplest one-dimensional case d = 1 (i.e. the answers to questions (B) 
and (C) are qualitatively different for these problems). Finally, we conclude the 
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paper by some remarks and open questions supported by numerical evidence, as 
well as the description of a similar problem we feel important for applications. 

2. The long-term problem 

The asymptotic behavior of the long-term optimal location Problem 1 1.1 1 has been 
intensively studied both using geometric (see e.g. [5] and [10] ) and variational meth- 
ods [3l 111], in the latter case mainly by means of T— convergence tools (see [8] for 
details on the theory). We summarize here the most important properties of this 
problem, in order to make later a comparison with the respective properties of the 
short-term problem. 

Notice first that 

(2.1) l n < Cn- 1/d . 

This estimate is straightforward, if one considers a set E composed by n points 
placed on a uniform grid of size approximately equal to n~ x l d . 

We use T— convergence theory to find the answers to the questions (A)-(C) posed 
in the Introduction. In order to apply it, we need to work with functionals defined on 
a common space, which we choose to be the space of all Borel probability measures 
D 5 (f2). To this aim we identify each set E C f2 having #E < +oo with the measure 
Pt. € CP(S1) defined by 

We define now a sequence of functionals on the space !P(£1) by setting 

V/ d F(E), if /x = Ms, #E < n, 



(2.2) £ n ( M ) := 



otherwise. 



The coefficient n x l d in the above formula prevents the minimization from degen- 
erating and is chosen according to (|2. 1[) . It is straightforward to recognize that 
minimizing $ n is equivalent to solving Problem 11.11 up to the above identification 
of sets with probability measures. 

We give here the T— convergence result only for the case v = f ■ £ . This result 
could easily be generalized to a generic measure v by inserting in the T— limit only 
the absolutely continuous part of v with respect to L d . The proof of this theorem, 
for the case when / is a lower semicontinuous function, can be found in 3 J . For the 
proof in the general case we provide below (i.e. when / S L 1 (f2) is not necessarily 
lower semicontinuous), one can apply the technique developed in [11] for the study 
of a similar problem (the so-called irrigation problem). 

Theorem 2.1. The sequence of functionals {5„} n Y -converges with respect to the 
weak* convergence of measures to the functional Soo- CP(O) — > K defined by the 
formula 



doo (m) : = 9d I -4/d d ^ d > where p = p ■ L d + p 



sing 1 



In P 

where p € L 1 (f2), p S ing stands for the singular part of the measure p, with respect to 
£j d , and 64 is a constant depending only on the dimension d satisfying < 8d < 00 
and given by 

9 d := inf \ liminf n l/d / dist (x, E„)dx : E„ C [0, l] d , #E„ < n > . 
{ n J[o.i] d J 

The most important consequences of this theorem are summarized in the follow- 
ing corollary and answer the questions we are interested in. 

Corollary 2.2. The following assertions hold. 
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(A),(B) One has 

Iimn 1 /^„ = min{^ 00 (/i) : [i G P(n)} = 9 d \\f\\ d/{d+1) > 0, 

hence, in 'particular, l n ~ n _1 / d . 
(C) Denoting by fi n :— /is„ ifte measures associated to a sequence of minimizers 
for Problem we have that fi n — 1 p, as v — > oo in the weak* sense 
of measures, where ft is the unique minimizer of 3m and is given by the 
formula 

n = cf d/{d+i). L d withc:= (J fd/(d+D dL d^ \ 

In particular, if v has constant density, i.e. v = c ■ L d , then jl has constant 
density as well, namely, ft = v . 

Proof. This statement comes from well-known properties of T— convergence (i.e. 
convergence of minima and of minimizers), once we find the unique minimizer of 
3oo- Due to the strictly decreasing nature of Joo with respect to p on {/ > 0} and 
to the fact that fJ, s i ng and pl{/ = o} do not affect the value of S^, it is straightforward 
that the minimizers should be absolutely continuous and concentrated on {/ > 0}. 
To identify the density of the absolutely continuous part, set A := f d /( d +i) . £<i L Q ) 
and w := pf~ d /{d+i) _ With this notation, finding the minimizers to ^Joo is equivalent 
to minimizing the functional w i— > J n w~ x / d dX over the set 

w G L 1 (fi, A), w > 0, / wd\ = l 



Jn 

By convexity of the map w i— > w^ 1 ^ d , it immediately follows from Jensen inequality 
that the minimum of the latter functional is attained at a constant function w. This 
shows fb = cf d /( d+1 ^ ■ L d and the computation of the constant c follows from the 
constraint fi G 3 3 (f2). At last, the value of min^oo is obtained by plugging /2 into 
the expression for □ 

The exact values of the constants 9d are known in the one-dimensional and two- 
dimensional cases. Namely, if 17 = [0, 1] and v = it is actually easy to compute 
explicitly the unique minimizer of Problem 1 1 . II which is given by the set of n points 
located at the centers of n equal disjoint intervals forming a partition of Q. In other 
words, one has 

^Ul^ 1 }' = sothat 01 = T 

i=i J 

In the two-dimensional case when il — [0, l] 2 and v = £ 2 lU it is known that the 
configuration of n points placed in centers of regular hexagons, is asymptotically 
optimal as n — > oo (see [S] or QH]), which gives possibility to compute explicitly the 
constant 02- Namely one gets 

|*l&= 31 y + 4 ^ 0.377 
1 ' 6\/2 3 3 / 4 

where a C R 2 stands for the regular hexagon of unit area centered at the origin. 

3. The short-term problem 

In this section we answer question (A) posed in the Introduction regarding the 
short-term optimal location Problem 1 1.21 Namely, we will show that 

(3.1) s n ~ n- x / d 

whenever v -C £ d , similarly to the asymptotic estimate l n ~ n~ x l d proved in 
Section H 
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Before proving (|3.ip we find it important to remark that the values s n actually 
may depend not only on n, but also on the chosen sequence of solutions {E^} of 
the short-term Problem [L2j This is due to the fact that at each minimization step 
both the position of the next point and the new minimum value may depend on 
the history i.e. on the configuration chosen on the previous steps. In other words, 
the choice of the optimal set at each step may affect the minimal values of all the 
following steps as the following example shows. 

Example 3.1. Consider the one-dimensional situation d = 1 with Q := [0,4] and 

v : = (2 ' 1[0 ' 1] + I 1 ! 2 ' 4 )) ' £l ' 

Then one clearly has that both singletons {1} and {2} are solutions to the short- 
term location problem at the first step n = 1, and both give the value s% = 5/4. 
Now, if one takes E^ := {1}, then at the second optimization step n = 2 we get the 
unique minimizer Yl 2 '■= {1, 3}, which gives the value S2 — 1/2. On the other hand, 
if at the first step one takes E^ := {2}, then at the second step one gets the unique 
minimizer Y,' 2 := {1/2, 2} which gives a different value s 2 = 5/8. 

Coming back to proving (|3.ip . we observe that it is impossible to use the T- 
convergence theory for this purpose. Namely, the constraint E^ +1 D ~E' n which is 
imposed at each minimization step, actually gives raise to a sequence of problems 
that, once rescaled and expressed in terms of probability measures as in Section [2j 
involve the functionals $' n given by 

{dn(fj), if either C = v^-i + 1 e 

or ^ = Me;_ 1 , 
+00, otherwise. 

It is therefore not difficult to see that whenever /Us' — 1 /2 in weak* sense of measures 
as n — > 00, then the T— limit functional would be finite only on /2 itself. This 
means that taking the limit of a sequence fj,s n and using the fact that it minimizes 
5'oo gives no additional information on the limit itself. Hence it is not possible to 
find p, in this way, contrary to the long-term case. We will therefore analyze the 
short-term Problem 11.21 directly, without using T— convergence tools. 
The following assertion is valid. 

Theorem 3.2. Let v < Z d . Then 

C in -^ d < s n < C 2 n- X l\ 

where the positive constants C\ and C2 do not depend on the choice of the sequence 
of solutions to the short-term optimal location Problem In particular, (|3.1|) 

holds. 

Proof. Let {T,' n } (resp. {S n }) be a sequence of minimizers for the short-term optimal 
location Problem 11.21 (resp. long-term optimal location Problem II. ip . where v = 
f ■ L d , f G i 1 (f2). Of course, the long-term cost is lower than the short-term one, 
namely, 

F(E„) < F(K)- 

This provides the required estimate from below as a consequence of Corollary 12.21 
It remains to prove an estimate of the form 

F(E' n ) < Bn- X l d 

for a suitable constant B independent of the choice of the sequence of solutions to 
the short-term optimal location Problem II. 21 

Once we have a set we want to estimate how much the functional F decreases 
when we add a point xq € fi. For the sake of brevity denote 5(x) := dist [x, E^). If 
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we set S := EJ, U {xo}, it is clear that for x £ B(xq, 5(xo)/4) we have dist (x, S) < 
J(xo)/4, <5(x) > j5(xo) and hence dist (x, S) < S(x) — S(xq)/2. Thus, if we set 

<5(x )\\ S(x Q ) 



g(x ) := v \B ^x , 
we get 

F(E) < F(S;) - . 9 (x ) and F(X' n+1 ) < F(V n ) - sup g(x ). 
To estimate sup Xo(EO g(xo) we use the inequalities 

g{x )dx = I dx [ dv{x)^-^-l{\ x _ Xo \ <s ^ Xo) / A }{x) 



n 



> 



> 



dv{x) [ ^^-^{\x-x \<&(x )/A}{xa)dx Q 
Jn 1 

d v ( x ) [ 5 S^Li 

{\x—x \<S(x)/5} 

(x )dx 

Jn 1 



d+l 



where tod stands for the volume of the unit ball in K d . In the above chain of 
inequalities we have used that, for xo S B(x, 5(x)/5) the condition |x — xo| < 
<5(xo)/4 is always satisfied as well as the inequality 5(xq) > §<5(x); moreover, the 
last inequality in the above chain is an application of Jensen inequality since v is a 
probability measure. 

We have therefore obtained 

sup g(x ) > i / g(x )dx >c( f dv(x)5(x)) = CF(Y>' n ) d+ \ 
x en Jn \Jn J 

which implies 

F(K +1 )<F(K)-CF(K) d+ \ 
where C > depends only on the dimension d of the underlying space. The conclu- 
sion follows now from Lemma 13.31 below (minding that the value Si is independent 
of the choice of the sequence of solutions to the short-term optimal location Prob- 
lemO). □ 



Lemma 3.3. Let {a n } be a sequence of nonnegative numbers satisfying a n +i < 
a n — Ca d+1 for all n £ N, where C > 0. Then there exists a number B > 
(depending only on a\, C and d), such that a n < Bn~ x / d for all n £ N. 

Proof. The proof will be performed by induction, simultaneously with the choice of 
B. 

The step n = 1 is satisfied choosing B > a\. Now we look for a condition 
on B such that the following statement is satisfied: if for some n £ N one has 
a n < Bn- 1 ^, then a n+1 < B(n + I)- 1 /*. Since 

a„+i < a n - Ca d n +1 < B n - 1/d - Ca d + l , 

then we must impose that 

Bn- 1/d - Ca d+1 < B(n + iy 1/d . 

Using again the fact that a n < Bn~ x l d , we get 

Bn-^ d < B{n + I)- 1 / 4 + CB d+1 n -^ d+1 ^ d , 

or equivalently 

(i + 

l/n 



' -< C B d 1 + - 
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Minding that 

(1 + ^-1 ,d(^) 



1/d, 



1/n dx 

it is sufficient to impose 1/d < CB d (f + l/n) 1 ^ , or even 1/d < CB d . Therefore, 
the choice B := ai V (l/dC) 1 ^ suffices for the proof to be concluded by induction. 

□ 

4. The one-dimensional case with uniform measure 

In this section we address questions (B) and (C) posed in the Introduction re- 
garding the short-term optimal location Problem 1 1. 21 In particular, studying ques- 
tion (B), we will show that, even in a very simple one-dimensional situation with v 
being the uniform measure on an interval, the ratio s n /l n (which is never smaller 
than one) does not tend to a limit as n — > oo, and hence, neither does n l ^s n 
(although by Theorem 13.21 one has n 1/d s n ~ 1), since by Corollary E2]lim„n 1 / d Z„ 
exists. Moreover, we will show that in this case even liminf„ s n /l n > 1. We further 
show that in the same situation there are infinitely many limit measures (in the 
weak* sense) of the sequences of /zjy as n — * oo, where Yf n solves Problem ll.2[ and 
what is more, neither of such limit measures is equal to the unique limit measure 
of the sequence of solutions to the long-term Problcm ll.il 

In this section we restrict ourselves to the case — [0, 1] and v — L 1 . We will 
further identify each set E C £1 having finite number of points with the partition 
of fl by the points of E. Namely, if E = {x\, . . . ,Xk}, we will always order the 
elements of E in such a way that 

X\ < X 2 < ■ ■ ■ < X k , 

and identify E with the partition of f2 given by the intervals (Aj)^ 1 , where 

Ai := [0,2:1], A 2 := [xi,x 2 ], ...,A k = [x k -i,x k ], A k +i ■= [x k , 1]. 

The intervals Ai and Afc + i will be further called external intervals, while all the 
other intervals of this partition will be called internal intervals. We will also identify 
the same set with the (k + l)-dimensional vector whose entries are the lengths of 
the intervals of the respective partition 

(|Ai|, . . . , |A fc+ i|) = (xi,x 2 - Xi,x 3 -x 2 ,...,x k - Xk-i, 1 - Xk), 

and we write 

S ~ (£1,2:2 - 2:1,2:3 -x 2 ,...,Xk- Xk-i,l - x k ). 

For instance, the set {1/3,2/3} is identified with the respective partition of [0, 1], 
and with the vector (1/3,1/3,1/3), i.e. {1/3,2/3} ~ (1/3,1/3,1/3). The intervals 
[0, 1/3] and [2/3, 1] are external intervals while the interval [1/3, 2/3] is internal. 

Proposition 4.1. Given a C f2, assume that E^ ~ (Ai, . . . , Afc+i), and let £fc+i 
be a minimizer of F over all sets E C Q such that E D Efc and ^E^+i = k + 1. 
Then E^+i = E^ U {2:}, while either of the two conditions hold. 

(1) Either x is the center of some internal interval Ai, so that 

Efe+i — (Ai, . . . , Aj_i, Xi/2, Xi/2, Xi+i, . . . , Afc+i), 

while 

F(Z k +i) = F(Eit) ~ A 4 2 /8, 

(2) or x divides some external interval Ai or A^+i with the length ratio 1 : 2 
closer to the boundary of [0, 1], so that 

Efe+i — (Ai/3, 2Ai/3, A2, . . . , Afe + i) or E^+i ~ (Ai . . . , Afe, 2Afc + i/3, A^+i/3) 

while 

F(E fe+1 ) = F(E fc ) - A?/3, 
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where i — 1 (if x *E Ai) or i = k + 1 (if x £ Afc +1 J 
Proof. If a; G Aj, i g" {1, fc + 1} (i.e. Aj is internal), then 

F{T, k+1 ) 



/ dist (z, £&) dz + / dist (z, {xi~i, x, Xi}) dz 

J[0,l]\Ai JA, 

/^ist (z^dz-^^"^ 
Jo 



2 

and to conclude, it is enough to note that the function 

[Xi-x,Xi] 3lH(lj-l)(l - £i-l)/2 

attains its maximum value A? / 8 at x = (xi-i + Xi)/2. The case when Ai is external, 
say, i — 1 (the case i = k + 1 is completely symmetric), is absolutely analogous, 
once one notes that 

F(S/ c+ i) = / dist (z, Sfc) dz + / dist (z, {x, xi}) dz 

J[0,l]\Ai Ja, 

" 1 x\ ( x 2 (x\ — x) 2 ' 



dist (z, E fe ) dz - } s } 
^ fc )-^ 2 " 2 ^ 



4 

Indeed, one has that the function [0,:ci] 3 x (3x 2 — 2x^x + x\)/A attains its 
maximum value A 2 /3 at x — a;i/3. □ 

We will further say that the point x G Ai is in optimal position, if either Ai is 
an internal interval and x is its center, or Ai is an external interval, and x divides 
it with the length ratio 1 : 2 closer to the boundary of [0, 1]. The above proposition 
says that whenever Efe c f2, #£fc = k and S^+i C fl, ^fcEfc+i = k + 1 solve 
Problem 1 1.21 then S^+i — U {a;} with a; in optimal position. 

We set now 



a := 



1 1 



2 • 3 i+1 ' 2 • 3* 



U 



1 1 
1 - -—,1 



2 • 3 4 ' 2 • 3 i+1 



i G N, 



so that, clearly gives a partition of f2, while |Qj| = 2/3 I+1 . This allows us to 

formulate the following corollary to the above Proposition 14. 11 

Corollary 4.2. Let S' fc C f2, #S' fe = k be a solution to Problem [LE Then, for the 
corresponding partition of f2 one has that for each interval Ai there exists a unique 
couple of numbers (j,h) G N 2 , such that |Aj| = '&~i2~ h , while 

(1) if Ai is internal, then j > 1 and Ai C f2j. 

(2) i/ Aj is external (i.e. i = 1 or i = fc + 1 i/ien h = 1 and 

f 0, fc = l, 

\ Jo or Jo ~ 1) otherwise, 

where jo := sup-jj : f2j D A TO /or some m = 2, . . . , fc}. 
Proof. The proof is easily obtained by induction on fc. □ 

Consider an arbitrary E' fe ~ (Ai, . . . , A^+i) solving Problem 11.21 Thanks to the 
corollary [472] we will identify each number A; corresponding to the internal interval 
Ai with the respective couple (j, h) G N 2 , j > 1, such that Ai = 3~i2~ h . We will 
further identify every external interval Aj with the couple (J + 1/2,-1/2), where 
j G N is such that Aj = 3 _:, 2 _1 (the reason for the latter identification will be 
explained in a moment). 

Consider now the set D of couples (j,h) G R 2 , where (j,h) G N 2 , j > 1, or 
j G N + 1/2, h = —1/2. We have that each interval A,;, and each Ai is identified 
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with a unique point d Xi S D. We introduce the ordering on D according to the 
following definition. 

Definition 4.3. We denote (j'2,/12) >~ (ji,hi), if 

h 2 -h 1 + (j 2 - ji) log 3/ log 2 < 0. 

Observe that the above relation well orders the set D. We have now the following 
simple assertion. 

Proposition 4.4. Let Si C ft, #S' fc = k be a solution to Problem \1.2l Letp,q 6 
{1, ...,k + l}, and consider d Xp = (j Xp , h Xp ), d Xq = (j Xq ,h Xq ) e D. Consider 
:= Yl' k U {x p } and £^+1 := T,' k U {x q }, where x p € A p and x q <E A q are in 
optimal positions. If d Xq >- d Xp , then 

Proof. One considers separately three cases. 

Case 1. Both A p and A q are internal intervals. Then 

- A 2 /8 

= F(£' k ) - 3- 2: >^2- 2hx p/8 

> F(Z' k ) - 3~ 2 ^2- 2h ^/8 = F(YT k+1 ), 

because h Xq -h Xp + {j Xq -j Xp ) log 3/ log2 < implies 3~ 2 ^p2- 2hx p < 3~ 2j ^ 2- 2hx * . 

Case 2. One of the intervals (say, A p ) is internal, another one (A g ) is external. 
Then, minding that h Xp = -1/2, and hence X p = 3~ jx P +1/2 2~ 1 = 3^ n p^ 1/2 2 2hx p , 
we get 

F(Zl +1 )=F(V k )-\ 2 p /3 

= F(T,' k )-3- 2j ^2 4h ^ 
= F(E' k )-3- 2j >>p2 ih *p 
>F(4)-3- 2 ^2- 2 ^/8 = F(E^. +1 ), 

because h Xt -h Xp +(jx q -j\ p ) log 3/log2 < 0andh Ap = -1/2 implies 3~ 2jx p 2~ hx p < 
3 -2ix < , 2 - 2ft ^/8. 

Case 3. Both intervals A p and A q are external, so that \ p = 3 _j *p +1 ' 2 2 _1 and 
A, = 3-^ +1 ' 2 2-\ while JXq < JXp . Then 

^ +1 )=F(£' fe )-A 2 p /3 

= F(E' fe )-3- 2 ^ +1 2- 2 /3 

> F(jy k ) - 3- 2 ^ +1 2- 2 /3 = Fm +1 ), □ 

As an immediate consequence of Proposition 14.41 we get the following theorem. 

Theorem 4.5. Let Y,' k C fl, = k be a solution to Problem \l. l A Let d = 

(j, h) G _D be maximal ( with respect to the order >- ) among the elements d Xq G D 
corresponding to q G {1, . . . , IfH' k+l C f2, #£J. +1 = fe+1 solves Problem \1.2\ 

then = U {x} where x € A p is in optimal position and p £ {1, . . . ,k + 1} 

is suc/i i/iai d = dA p (such an element may be nonunique) . 

We are able now to formulate the following statement which says exactly how 
the sequence of solutions to Problem 11.21 looks like in the particular case we 

are considering. 

Corollary 4.6. Enumerate D in the order given by the relation y, so that D = 
{d(i)}iZi, d{i + 1) y d(i). The sequence of solutions £ fe C fl, ff£' k = k to 
Problem can be described in the following way by induction on (D,y). Each 
d{i) = h(i)) £ D gives rise to a part § l of the sequence of optimal sets. 
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• §° consists of the unique set ~ (1/2; 1/2). 

• The set S l+1 consists of the sets , l5 . . . , where S' fe is i/ie Zasi element 
of&, 

S i+i = Sj-U{xj}, j = k,...,h- 1, 

h := k + 2' 1 W — 1, .Tj G A p C fij(i) in optimal position, and p G {1, . . . , j} 
is an arbitrary index satisfying d{i) = d\ p . 

Proof. The proof is easily obtained by induction on the well-ordered set (D,>~). □ 



Remark 4.7. An easy consequence of what has been proven so far, is that in the 
case we are considering (that is, when v is the uniform measure over the interval) 
the value of s n does not depend on the particular sequence of solutions {£^} to 
short-term Problem 11.21 i.e. it depends only on the index n. 

Example 4.8. We consider the first elements of the possible sequences of solutions 
E' k C ft, #E' k = k to Problem O 

(1) i = 0, S° = {S'J, where 

Ei~ (1/2,1/2). 

(2) i = 1, d(l) := - (-1/2,-1/2) G £>. From d(l) = d Ap we get 
A p = 3 _ J'( 1 )+ 1 / 2 2 _1 = 1/2, which corresponds to the two external intervals 
of length 1/2. Hence, S 1 := {£' 2 ,£^}, where 

either £' 2 ~ (1/6, 1/3, 1/2) or £ 2 ~ (1/2, 1/3, 1/6), 
(1/6,1/3,1/3,1/6). 

(3) i - 2, d(2) := (j(2),/i(2)) = (1,0) G D. From d(2) = d Xp we get A p = 
3-i(2)2~' l ( 2 ) — 1/3, which corresponds to the two internal intervals of length 
1/3, both belonging to ftj(2) = ^i- Hence, S 2 := {£ 4 , Eg}, where 

cither £4 ~ (1/6, 1/6, 1/6, 1/3, 1/6) or £ 4 ~ (1/6, 1/3, 1/6, 1/6, 1/6), 
E' B ~ (1/6, 1/6, 1/6, 1/6, 1/6, 1/6). 

(4) i = 3, d(3) := (j(3),/i(3)) = (1/2,-1/2) G D. From d(3) = d Ap we get 
A p = 3 _J ' 3 ' 1 / 2 2 _1 = 1/6, which corresponds to the two external intervals 
of length 1/6. Hence, S 3 := {E' 6 ,E' 7 }, where 

either £ 4 ~ (1/18, 1/9, 1/6, 1/6, 1/6, 1/6, 1/6) 
or £ 4 ~ (1/6, 1/6, 1/6, 1/6, 1/6, 1/9, 1/18), 
£' 5 ~ (1/18, 1/9, 1/6, 1/6, 1/6, 1/6, 1/9, 1/18). 

(5) i = 4, d(4) := (j(4),/i(4)) = (1,1) G £>. From d(4) = d Ap we get X p = 
3-i( 4 )2 _ ' l ( 4 ) — 1/6, which corresponds to the four internal intervals of length 
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1/6, all belonging to £ljU) = ^i- Hence, S 3 := {E 8 , Eg, E' 10 , S^}, where 
either E 8 ~ (1/18, 1/9, 1/12, 1/12, 1/6, 1/6, 1/6, 1/9, 1/18) 
or E 8 ~ (1/18, 1/9, 1/6, 1/12, 1/12, 1/6, 1/6, 1/9, 1/18) 
or E 8 ~ (1/18, 1/9, 1/6, 1/6, 1/12, 1/12, 1/6, 1/9, 1/18) 
or E 8 ~ (1/18, 1/9, 1/6, 1/6, 1/6, 1/12, 1/12, 1/9, 1/18), 
either E 9 ~ (1/18, 1/9, 1/12, 1/12, 1/12, 1/12, 1/6, 1/6, 1/9, 1/18) 
or E 9 ~ (1/18, 1/9, 1/12, 1/12, 1/6, 1/12, 1/12, 1/6, 1/9, 1/18) 
or E 9 ~ (1/18, 1/9, 1/12, 1/12, 1/6, 1/6, 1/12, 1/12, 1/9, 1/18) 
or E 9 ~ (1/18, 1/9, 1/6, 1/12, 1/12, 1/12, 1/12, 1/6, 1/9, 1/18) 
or Eg ~ (1/18, 1/9, 1/6, 1/12, 1/12, 1/6, 1/12, 1/12, 1/9, 1/18) 
or Eg ~ (1/18, 1/9, 1/6, 1/6, 1/12, 1/12, 1/12, 1/12, 1/9, 1/18), 



We are able to claim now the following assertions answering questions (B) and (C) 
for the case we are considering. 

Theorem 4.9. Assume v — £ 1 l[0, 1]. Then the bounded sequence {ns n } does not 
converge to a limit as n — > oo. Further, for any sequence {E^} of solutions to 
Problem the sequence of probability measures £isy has infinitely many limit 
measures in the weak* sense as n — > oo. 

Proof. Fixed a j £ N, consider the sequence of indices k n such that EJ, induces 
a subdivision of the set flj into 2™ equal subintervals while the partition of £1 
corresponding to E' fe +1 divides further one of the latter subintervals into two equal 
intervals. This means that the sequence of partitions corresponding to optimal sets 
Ti' k , starting from k — k n and up to k = k n + 2™, will be obtained by dividing into 
two equal parts at each step one of the 2" subintervals in Qj. Set now 

a n := / dist (x, E^ )dx and b n := / dist (x, E^ )dx, 

./[0,l]\fij Jttj 

Fix an arbitrary dyadic number A £ (0, 1) and consider the new sequence of indices 
k' n := k n + A2™, so that the partition of corresponding to E^ coincides with that 
corresponding to E/^ up to the fact that some of the subintervals in fio, namely, a 
fraction A of the total, have been split in two parts. This implies that 

k'J K = k' n F{T,' K ) = (k n + A2")(a„ + (1 - A/2)6„), 

splitting an interval in two equal parts reduces the average distance from its points 
to E by a factor two. Mind that 

k n lk„ = k n F(Ti kn ) — k n {a n + 6 n ), 

and assume by contradiction that, lim^ klk exists. Then lim„ k' n lk' n = \\va n k n lk nl 
and hence, 

{kn + A2")(a„+(l - ^)fo n ) - k n (a n + b n ) 

= \(2 n a n -^k n b n + fl-^) 2 " fe ™ 
X(2 n a n ~h n b n + l(l-^)) ^0 



2 9 V 2. 

as n — > 00. Therefore, 2 n a n — k n b n /2 — * A/18 — 1/9 as n — > 00, which means that, 
taking two different dyadic values of A £ (0, 1), the sequence {2 n a n — k n b n /2} has 
two different limits. This contradiction proves the first claim. 
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To prove the second claim, we first show that the limit measures of the sequence 
{/isj} are not unique, where / € [k n , k n + 2 n ]. To this aim, suppose by contradiction 
that all the subsequences of the above sequence converge in the weak* sense to the 
same limit measure /i as n — > oo. In particular, this is the case of {a 1 ^ }, which 
clearly then converges to a uniform measure over Qj since the points of Y,' k arc 
uniformly distributed over Qj. Observing that k n > 2™, and hence 

1/(1 + A) < k n /k' n < 1, 

we may assume without loss of generality that up to extracting a subsequence of 
k n (not relabeled), the sequence {k n /k' n } converges to some finite limit. Hence 

1 = lim ^7T\ = lim = (! + A ) lim TF> 

n— »oo fj,^ ("jj n— >oc Z n / K n n—>oc ft' n 

which gives lim n k n /k' n ^ 1. But since 

then passing to a limit in the above relationship as n — > oo, we get 

k X2 n 
fi(Qj) = lim -rj-K^j) + lim -rr 

n n 

= n(flj) lim lim (l-^ 

n^oo ft' n^oo y fc„ 

Minding that lim„ k n /k' n ^ 1, the above relationship is only possible when /x(f2j) = 
1, i.e. when the unique limit measure [i is concentrated on Qj, which is clearly not 
the case because 



(4.1) //(% +1 ) > 



fi(Qj) 



3 

To verify the latter inequality, it is enough to notice that 

(%+i) 2 h + 1 . , . . . . . . 

Ms' fer> ("j) 2™ + l 

so that h > n — log 3/ log 2, which gives (|4.ip in the limit as n — » oo. 

At last, to show that the limit measures of the sequence {/^s^} are infinite, it is 
enough to vary j G N. □ 

In Figure [1] below we provide a graph of s n /l n for the case we are considering 
calculated according to the algorithm given by Theorem 14.51 and Corollary! 



Theorem 4.10. No sequence of solutions to Problem (with v uniform 

measure on [0,1],) has the Lebesgue measure rC^fO, 1] as a limit measure of some 
subsequence of /is' in the weak* sense as n — > oo. 

Proof. By using Corollary 14.61 to any k we can associate a set Qj in the following 
way. Let i £ N be such that T*' k E S\ Consider now the pair (j(i),h(i)) € Z? 
corresponding to the index i, and set j :— j(i). The index j = j[k) associated to a 
set Ti' k represents the set Qj where we last put a point in building and also the 
set Qj which is being divided in a dyadic way by the optimal sequence at step k. 
Namely, if we look at the points composing S' fc , we see the following picture. 

• For every j > maxk'<k j[k'] no point of belongs to Qj. 

• Every Qj with j < maxt'<t j[k'] and j ^ j[k] contains a certain number 2 h 
of points of T,' k (precisely we have h = max{/i : (j, h) -< (j[k], k)}). 

• At last, in Qj[k] we cannot exactly predict the number of points of EJ.. 
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Figure 1. Plot of s n /l n for the case of uniform density over [0, 1]. 

Up to choosing a subsequence of n (not relabeled), we can suppose the existence of 
two distinct indices j\ ancl 3% such that for any n (from the chosen subsequence) 
we have j[n] ^ j\, and j[n] ^ ji. Extracting a further subsequence of n (again not 
relabeled), we may assume that maxfe'<„ j[k'\ > (j\ V j'2), and hence the number of 
points of T/ n in the sets and flj 2 is a power of two. This means that, if we set 
Mn := ^T,' k and we suppose fi n — fi = -C 1 l[0, 1] in the weak* sense as n — > 00, we 
get 

2 A(n,ii,j2) _ _^ _ gis-ji 

Mn(%) M(%) 

where A(n, 31,32) is an integer exponent depending only on n, j\ and j2- Yet this 
is a contradiction as a sequence of powers of two can converge only to 0, +00 or a 
power of two. □ 

The above theorem gives a rigorous formulation of the insight of the authors 
of [15] who noticed that a short-term strategy in the above one-dimensional sit- 
uation, compared to other allocation policies, induces "a more uniformly spaced 
allocation yet not completely unijorm" . 

Corollary 4.11. One has 

1 < liminf s n /l n < limsups n /Z n , 

n >oo n — » qq 

if v is the uniform measure over [0, 1]. 

Proof. The second inequality is just the reformulation of the first claim of Theo- 
rem 2jJ] minding Corollary 12.21 To prove the first inequality, note that s n > l n) 
while, should the above liminf be equal to one, we could build a subsequence 
{Sjj } of solutions to short-term Problem 11.21 such that the respective sequence of 
measures {ms^ } is asymptotically optimal for the sequence of functionals {$k n }n 
defined by (12. 2p . From general results in T— convergence theory we would obtain 
/i£ fcm — v fi in the weak* sense as n — > 00, where p, is the limit measure for the 
long-term Problem 11.11 i.e. the Lebesgue measure over [0,1], which contradicts 
Theorem OHl □ 

5. Concluding remarks and open problems 



We conclude the paper by a list of remarks and open problems regarding the 
model studied in this paper, which we consider to be interesting for further study. 
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Question I. It seems interesting to understand whether the phenomena stud- 
ied in Section 0] (e.g. non existence of the limit of s n /l n , the fact that the re- 
spective liminf is strictly greater than one, non- uniqueness of the limit measures 
for the short-term location problem etc.) for the case v = £ 1 l[0, 1] occur in the 
case of generic measure v in one-dimensional case d = 1. Further, it seems to 
be important for applications to obtain sharp estimates on liminf n^oo s n /l n and 
limsup n _ >o0 s n /l n even for the uniform measure v (in the latter case, the numerical 
computation provided in this paper (Figure [1]) suggests that both limits are close 
to 1, so that, in practice, one could assume that asymptotically the value of the 
functional in the short-term problem is almost the same as in the long-term one, 
which would give a rigorous statement of the idea first suggested in [TS]) and to 
characterize completely the limit measures of the short-term problem. 

Question II. In the case of the generic space dimension d > 1 we are able to 
prove only the asymptotic order estimate of the value of the functional in the short- 
term problem. It seems therefore important to study the above posed problems 
for generic space dimension even in the simplest case when, say, v is the uniform 
measure on the unit ball or on the unit square. Figure [5] shows what one can expect 
about the behavior of s n in the case of a uniform density on a unit square (for conve- 
nience, the normalized values s n /l^ are provided, where := n~ 1 / d 62\\f\\d/(d+i) 
stands for the asymptotical value of the long-term minima according to Theo- 
rem 12.11) , while Figure [3] shows how the distribution of points in this case looks 
like. Both figures are obtained by a numerical calculation on a uniform 200 x 200 
grid. More results of numerical computations for the short-term optimal location 
problem can be fund on the web page [18j . 
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FIGURE 2. Plot of s n /Z£° for the case of a uniform density on a unit square 

Question III. At last, it is natural to mention here a similar problem intro- 
duced in [5] and sometimes called irrigation problem (see [H |H1 I12[ 1131 114j ) on 
minimization of the average distance functional but over compact connected sets of 
finite length rather than over discrete sets of points like in this paper. The state- 
ment of such a problem is obtained by replacing the constraint on cardinality 
of the unknown minimizer £ in the location problem by the similar constraint on 
the one-dimensional Hausdorff measure 3-C 1 (S). The problem is therefore that of 
finding a minimizer of the cost f:EH J n dist (x, S) dv{x) over all compact and 
connected sets S C f2 satisfying 5{ 1 (E) < I with given I > 0. This problem, which 
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Figure 3. Short-term location for the case of a uniform density 
on a unit square of (left) the first 9 points (right) 1000 points. 

again can be interpreted as the long-term one, admits also a natural short-term 
approach. Namely, by letting the parameter I increase, we would like to find an irri- 
gation set £ which increases continuously with I, in such a way that the functional 
F decreases as fast as possible. This construction may be made rigorous through 
a slight modification of the well-known method of minimizing movements (see for 
instance [TJ [2] for the presentation of the theory) . Namely, fixed an arbitrary time 
step r > 0, we minimize the functional £ h- ► F(£) + 3i 1 (E \ ££) 2 /2r over the set 
of all connected compact subsets of £1 satisfying the constraint ££ +1 3 ££, where 
£J := (alternatively, one could drop the JC 1 — penalization term in the above func- 
tional and add the additional constraint JC 1 (£*,+! ) < Jf 1 (Ep + t instead). Let 
then 1 i— > £ r (i) be the piecewise constant map defined by 

57(f) := ££ for t G [kr, (k + l)r[. 

One should first study whether in this way one can find a well-defined increasing 
evolution £(i) as a limit (in the Hausdorff topology) of £ T (i) as r — * 0. Further, 
several questions then arise about the behavior of £(i). For instance, it has been 
proven in [6l El [T3l [14] that the solutions of the long-term irrigation problem under 
suitable conditions on problem data do not contain loops, have a finite number of 
endpoints and may have only regular tripods (i.e. triple junctions with the branches 
with infinitesimal angles 120° between each other) as branching points, which are 
at most finite in number. It is interesting to verify whether the same or similar 
properties hold also for solutions of the short-term problem E(t) for all (or for 
some) t € R + , or whether such solutions are just simple curves without branching 
points. Besides, one is also interested in the asymptotic density of the curve £(£) 
(i.e. the weak* limits of the measures , K 1 L.T,(t)/'K 1 ('S(t)) as t — > oo), which for the 
long-term irrigation problem was studied in [TT] . 
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